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Models of late-time neutrino mass generation contain new interactions of the 
cosmic background neutrinos with supernova relic neutrinos (SRNs). Exchange of 
an on-shell light scalar may lead to significant modification of the differential SRN 
flux observed at earth. We consider an Abelian U(l) model for generating neutrino 
masses at low scales, and show that there are cases for which the changes induced 
in the flux allow one to distinguish the Major ana or Dirac nature of neutrinos, as 
well as the type of neutrino mass hierarchy (normal or inverted or quasi-degenerate). 
In some region of parameter space the determination of the absolute values of the 
neutrino masses is also conceivable. Measurements of the presence of these effects 
may be possible at the next-generation water Cerenkov detectors enriched with 
Gadolinium, or a 100 kton liquid argon detector. 



INTRODUCTION 



Neutrino flavor conversion has been observed in the solar (SuperK, SNO) atmospheric 
(SuperK) [^J, and terrestrial (KamLand, K2K) 3] neutrino data, providing evidence for non- 
vanishing, sub-eV neutrino masses. There now remains the longstanding theoretical question 
of how the neutrinos acquire their masses. The most elegant solution to this puzzle is the 
seesaw mechanism (4]: one assumes that lepton number is violated at some high scale Al in 
the form of right-handed neutrino, N, Majorana masses, Mn ~ Al. This induces, at a lower 



2 



scale, an effective operator of the form 0(1) x (LH) 2 /A^ , where L denotes a lepton doublet 
and H the Higgs field. The oscillation data then imply that Al ~ 10 14 GeV . However, it is 
difficult to devise an experimental test of this mechanism (see however Q]). Therefore, it is 
important to explore alternate natural mechanisms for neutrino mass generation, especially 
those that may be tested in experiments at low energies. 

A class of such models that have astrophysical and cosmological tests are the models of 



late-time neutrino mass generation 



, 13, |8j . In these models, neutrino masses are protected 
by some flavor symmetry different from the one related to the charged fermion masses, for 
example some global £7(1) at symmetry. The small neutrino masses are generated when the 
new symmetry is broken at low scales. The effective Lagrangian for these models can be 
schematically written for either Dirac or Majorana particles where the neutrino fields are 
neutrino mass eigenstates as 

C° v = C km + y u <f)uN + V {(/)), C K J = C kin + y v <j>uu + V{$), 

where C^n is the kinetic piece of the lagrangian, y v is a dimensionless coupling, v is a 
standard model neutrino field, N is an extra field introduced for the case of neutrinos being 
Dirac particles, <fi is the scalar field and V(<f>) is the associated scalar potential. After 
spontaneous symmetry breaking the neutrinos acquire masses given by 

m v = yj (1) 

where m u is the mass of a particular neutrino mass eigenstate and / is the symmetry 
breaking scale (/ = (0), where (0) is the vacuum expectation value (VEV) of 0). With just 
one scalar the couplings are diagonal in the mass basis. In addition, since the scalar couples 
to neutrinos only, the constraints on the symmetry breaking scale / are weak 

In addition to generating neutrino mass through their VEVs, the new light scalars provide 
another neutrino-neutrino interaction process aside from the Standard Model Z° exchange. 
The effects of the neutrinos coupling to these scalars on the cosmic microwave background 
has been previously studied Constraints have been placed on the symmetry breaking 
scale, /, the scalar mass, Mq, and the scalar-neutrino couplings, y v , in these models from 
cosmological considerations [3, ls[ [9], as well as from demanding that supernova cooling 
and the flux of the 1987a neutrinos would not be significantly modified in the presence 
of the additional fields 
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111] (for constraints related to generating the observed baryon 



asymmetry of the universe see 
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In this paper we show that the presence of this new physics significantly modifies the 
spectrum of supernova relic neutrinos (SRNs) at earth. This modification occurs because 
SRNs can interact with cosmic background neutrinos through exchange of the new light 



scalar. In [Uj this effect was studied assuming a single flavor, Majorana, case. Here we 
extend our study and include various interesting aspects related to the nature of the neutrino 
flavor sector, for example the breaking of lepton number, the effect of multiple generation, 
etc. For simplicity we confine our study here to late neutrino mass models with a single 
£7(1) symmetry. We show that the energy spectrum of the SRN flux is sensitive to the 
type of neutrino mass hierarchy and whether neutrinos are Majorana or Dirac particles. We 
discuss how in some specific cases one can get additional information about the neutrino 
masses as well. In addition, detection of this signal would also be direct evidence of the 
presence of the cosmic background neutrinos. Hundreds of events per year from the flux 
of the SRN antineutrinos could be seen at next-generation large megaton water Cerenkov 



detectors 13l. Il4| such as UNO, Hyper-Kamiokande or MEMPHYS if they are enriched with 
Gadolinium jisj j. or from the flux of the SRN neutrinos in a large 100 kton liquid argon 



neutrino detector 



In section [TT] we show how the SRN flux, including cosmological evolution, is modified 
through the new interactions. In section II I II we consider the normal and inverted neutrino 
mass hierarchy cases as well as quasi-degenerate neutrino masses. We also consider the 
possibility of neutrinos being either Majorana or Dirac particles. In addition we show that 
there is a particularly interesting signal that leads to the determination of ratios of neutrino 
masses. The conditions for establishing a statistically significant signal above background 
are discussed in subsection G of section III. Finally, in section HVt we conclude. 

II. THE SUPERNOVA RELIC NEUTRINO FLUX 

The resonance interaction of the SRN with cosmic background neutrinos through the 



exchange of a new light scalar was previously discussed in ll| for a single neutrino mass 
eigenstate with m u ~ 0.05eV. 

In this paper we consider the case of the scalar interacting with three neutrino mass 
eigenstates and show the very interesting effect on the observed SRN spectrum for the 
normal mass hierarchy, inverted mass hierarchy, and for quasi-degenerate neutrino masses. 
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We briefly discuss the case of adding one sterile neutrino. 

We start with the SRN flux without the new interactions. The diffuse SRN flux is a 
remnant of neutrinos emitted from all the supernova that have occurred in the universe [ij 
This flux is given by 

dt 



z max dN({\ + z)E v ) 

j^sn{z) — 



1 + 2 



dz 



dz 



(2) 



10 dE u 

where Rsn{ z ) is the comoving rate of supernova formation, dN((l + z)E u )/dE u is the 
neutrino energy spectrum emitted by supernova, dt/dz is for the cosmological expansion, c 
is the speed of light, z is the redshift, and E v is the neutrino energy. 

The quantity R$n is the comoving rate of supernova formation, which can be parame- 
terized as [17] 

/n rm\ 

(3) 



Rsn{z) 



P* \z), 



where p* (z) is the star formation rate given by 



(4) 



p\ ( z ) = (1 - 2) x 10" 2 M Q yr^Mpc' 3 x (1 + zf. 

We take R SN (0) = 2 x lO^yr" 1 Mpc^ 3 , (3 = 2 (for < z < 1) and (3 = (for z > 1) [Tl 
These are 'median' values for the parameters which have uncertainties in them coming from 



the uncertainty in the present knowledge of the cosmic star formation rate 
factor dt/dz is given by 



13, 



dt 

dz 



100- 



km 



h (l + z)y/n M (l + z) 



IB. 1 The 



(5) 



s Mpc 

with n M = 0.3, Vl A = 0.7 and h = 0.7. 

The energy spectrum of the neutrinos emitted by a supernova has been modeled by several 



groups [19|, |20|, |2lj . One of the models leads to perfect equipartition of the energy radiated 
into each neutrino flavor [lj| . A second model makes a detailed one dimensional calculation 
of all relevant neutrino processes in the collapsing star and uses a variety of supernova 



progenitor masses [21|J. Another model proposed by Keil, Raffelt, and Janka [20] (KRJ 
model) performs their calculations using a MC simulation. These models give spectra which 
have a width narrower than that of a thermal spectrum (so-called "pinched" spectrum 22[). 



It also predicts an average energy for the muon and tau flavor neutrinos very close to 



1 Future SN observatories will have the power to significantly reduce the related uncertainties [14] 
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the average energy of the electron antineutrinos. As an illustrative example we take the 
neutrino spectrum given by the KRJ model with the additional assumption that the total 
energy carried by each neutrino flavor is L Ue = L Pe = L Ux = 5 x 10 52 ergs, where x stands 
for the muon and tau neutrinos and antineutrinos. The KRJ energy spectrum of neutrino 
flavor eigenstates produced by a supernova is given by 

dN u (E y ) (1 + B V ) 1+ P"°L V (E v V Va r , „ .E v . 

exp[-(l + #,J=]i (6) 



dE »~ r[i + p Ua ]E Ua 2 \E V J a E Ua 



with the E„ representing the average neutrino energies and (3 Ua characterizing the amount 



of spectral pinching. Here for simplicity we set the values of the parameters to 20j (for 



recent numerical studies see e.g. [231 ]) 





= 3.4 




= 13.0 MeV, 


Ve ■ fip e 


= 4.2 




= 15.4 MeV, 




= 2.5 




= 15.7 MeV, 



(7) 



and also neglect effects such as shock wave and turbulence 24|. 

Because of matter oscillation effects, neutrinos emerge from a supernova as coherent fluxes 
of mass eigenstates which we label as F Vi , where % — 1, 2, or 3 represents the particular 



neutrino mass eigenstate 



25). 



If neutrino flavor evolution inside of the collapsing star is either fully adiabatic or fully 
non-adiabatic (the flavor evolution is adiabatic if the mixing angle sin 2 8i3 > 10 -3 and non- 
adiabatic if sin 2 8i3 < 10 -5 ) then the energy spectrum of each neutrino mass eigenstate 
that leaves the surface of the star corresponds to the original energy spectrum of some 
particular neutrino flavor eigenstate at emission from the neutrinosphere, i.e., there is a 
one-to-one correspondence between each dN Ua (E Ua ) j dE Va and some dN Vi (E Vi ) / ' dE Vi ? The 
original produced flux of some neutrino flavor at the neutrinosphere will be labeled as F® a . 
Translated back into the flavor basis, the expressions for the v e and v e fluxes emerging from 
a supernova can be written as 

F Ve = P H \U e2 \ 2 Fl + (1 - P H \U e2 \ 2 )Fl , (8) 
F Pe = \U el \ 2 Fl + \U e2 \ 2 F° 



See, for example, Table 1, Fogli, et. al. 
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for the normal mass hierarchy and 



\U e 2\ 2 F® a + |t/" e i| 2 F° 



F Ve = P H \U el \ 2 F?+(l-P H \U el \ 



(9) 



for the inverted mass hierarchy, where P#(and Ph) = for the adiabatic case and 



Pf/(and Ph) = 1 for the non-adiabatic case 23]. In the equations above, \U, 

COS 2 6>i 2 , \U P o\ 2 



sin 2 # 12 , 6j2 



0.86 ±0.3, and \U e3 \ 2 « 



9 e (0 Q is the solar mixing angle 26J), where sin 2 IQyi 



>el\ 



271 ]. When supernova neutrino flavor evolution is non-adiabatic 



then the v e and P e flux for the normal and inverted hierarchies are identical. 

We show that the addition of a new light scalar opens the possibility of determining the 
neutrino mass hierarchy independent of the neutrino fluxes, and also independent of the 
adiabatic or non-adiabatic nature of supernova neutrino flavor evolution. 



A. Modifications due to New Physics 

We consider the modifications to the SRN flux due to the resonance interaction of an SRN 
with a neutrino in the cosmic neutrino background. In this process a supernova neutrino 
with energy E„ will go through the resonance when the kinematic condition 

is satisfied. More specifically, a neutrino observed with energy E„ will have gone through 
resonance if its energy lies in the region 

rpRes 

~7iObs ^ rpRes 



< E^ bs < E^ s , (11) 



1 + z 

where z is the redshift. There will be large depletion of the SRN flux in the energy domain 
given in Eq. [11] as long as the neutrino-scalar coupling satisfies 11] 



y > 46 x 10 ~ 8 T^v ■ (12) 

This condition comes from requiring that the mean free path for absorption is much smaller 
than the Hubble scale. It is important to note that in the narrow width approximation for 
the resonance, the condition Eq. [12] is a sufficient condition to guarantee the absorption of all 
three neutrino flavors. After the neutrinos that have energies in the region given by Eq. [11] 



go through the resonance they are redistributed to lower energies when the produced scalar 
decays back to neutrino mass eigenstates j^. In particular, neutrinos after interaction will 
be redistributed with a flat energy distribution from zero energy up to the original energy 
of the incident supernova neutrino. 

To find the effect of these interactions on the flux of the SRN we note that the neutrinos 
leaving a supernova at redshift z emerge as the mass eigenstates. However, these mass 
eigenstate fluxes are now modified through interaction with the cosmic background neutrinos 
as they propagate to the earth. We consider for simplicity an Abelian U(l) late neutrino 
mass model. This implies that the Yukawa interaction between the scalars (in particular 
the Goldstone) and the neutrinos are diagonal in the mass basis (this is not the case in a 
model with non- Abelian symmetries). Supernova neutrinos are in their mass eigenstates and 
each mass eigenstate interacts only with the same mass eigenstate background neutrino via 
Goldstone exchange. To illustrate how the interactions modify the neutrino mass eigenstate 
flux we consider as an example the flux of the v\ mass eigenstate, F Vl . 

We start by defining the modified flux of the V\ eigenstates as F Vl . For each redshift z 
the v\ eigenstates that satisfy the condition given by Eq. [TT]will have resonance interaction 
with cosmic background neutrinos, producing the intermediate scalar. A neutrino mass 
eigenstate will go through the resonance when the coupling satisfies Eq. [T2j The cross 
section (averaged over the width of the resonance) for this to occur is approximately given 
by ORes — k/Mq. This will lead to a mean free path much smaller than the typical distance 
a supernova neutrino will travel to arrive at the earth, for the values of Mq that we consider. 

We label the absorbed flux as F^ es . Naively, the modified flux would be given by 

F V1 = F U1 -F^ S . (13) 

However, this expression does not take into account that the scalar decays back into neutrino 
mass eigenstates. We need to add this contribution to Eq. [13j The scalar can decay to any of 
the neutrino mass eigenstates. The probability that the scalar decays to a particular neutrino 
mass eigenstate is proportional to the square of the Yukawa coupling of that particular 
neutrino mass eigenstate to the scalar. From Eq. [I] we note that the relative probabilities 
are proportional to the ratios of squares of the neutrino masses, 



s 



Then the probability that a scalar decays to the neutrino mass eigenstate v\ is P\. These 
decays result in redistribution of the neutrino energies from zero energy up to the energy of 
the incident SRN with a flat energy distribution. We define Pi x Ff_^\, as the fraction of 
the flux of v\ that initiate a resonance, producing a scalar which then decays back into a v\ 
eigenstate with degraded energy (indicated by the notation 1'). Then, Eq. [13] is modified to 

F~ = F U1 - F^ s + P,x F*? v . (15) 

We still need to take into account the contributions from the decays of scalars produced 
by other neutrino mass eigenstates. Therefore, there should be a sum over all of the initial 
states, and Eq. [TBI becomes 

F V1 =F V1 -F^ s + P l Yl F ^v- ( 16 ) 

i=l, 2,3,1,2,3 

In more general notation, for the j th neutrino mass eigenstate, 

t=l,2,3,l,2,3 

The contributions over a range of redshift must be taken to determine the total flux at 
earth. If neutrinos are Dirac particles then there is factor of 1/2 multiplying the last term 
(see discussion in Section UlI Ej) . 

The modified flux of electron neutrinos and electron antineutrinos can then be written as 

FZ e = cos 2 6 l2 F Ul + sin 2 6 l2 F U2 , (18) 

and 

% e = cos 2 OnK, + sin 2 e l2 % 2 . (19) 

Finally, we note that each neutrino mass eigenstate goes through resonance at different 
energies given by Eq. [10] when there is just a single scalar of mass Mq- Depending on the 
details of the neutrino mass hierarchy, these resonance energies can either be very close to 
one another, or widely spaced apart. 

III. SIGNALS OF MODELS OF LATE-TIME NEUTRINO MASS GENERATION 

In this section we discuss the signals for the neutrino mass hierarchy in the observed SRN 
flux. We consider the case of the normal neutrino mass hierarchy, the inverted neutrino mass 



E 



nijes 



(17) 
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hierarchy, and the possibility that the neutrino masses are quasi-degenerate. We also show 
the effects of the neutrinos being Dirac or Majorana particles on the SRN flux signal. For 
these cases, unless otherwise noted, we choose the value of E^ es , defined in Eq. [101 at 
z = to be equal to 15 MeV for one of the neutrino mass eigenstates. This choice is 
made to illustrate the effects of the resonance process and to determine a region of the 
parameter space of the late-time neutrino mass generation models where the effect of the 
SRN modification would be seen. This is so, since for water Cerenkov detectors near reactors 
the background becomes negligible above about 13 MeV Q, Q. As will be discussed in the 
conclusions, the energy resolution is about 2 MeV. In subsections A through E we consider 
the case where neutrino flavor evolution in the SN is adiabatic, and in subsection F we show 
that the same features are obtained for a case where the flavor evolution is nonadiabatic. In 
subsection G we discuss the detection of the new interactions. 

In the following we focus on the flux of electron antineutrinos that arrive at earth, since 
the proposed water Cerenkov experiments for detection are sensitive to this neutrino flavor 
through the interaction of electron antineutrinos with protons with a cross-section given 



by |29| 



a = 10-43 pe+jBe+jB; 0.07056 + 0.020181n^-0.0019531n^, cm 2 ; (2Q) 

where E e + and p e + are the energy and momentum of the detected positron. Note that 



detection of the electron neutrino component o 



would provide complementary information 



30]. 



the SRN flux at a large liquid argon detector 



A. Normal Neutrino Mass Hierarchy 

As an example to illustrate how resonance interactions between the SRN and the cosmic 
background neutrinos can affect the SRN flux, we first consider a normal mass hierarchy of 
neutrino mass eigenstates. As a particular example of this hierarchy we choose the masses 
of the mass eigenstates to be 

mi = 0.002 eV, 

m 2 = 0.009 eV, (21) 
m 3 = 0.05 eV. 
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'orms to the best value of the atmospheric mass splitting \m^ — m\ 2 \ — 2.4 x 



This con 

10~ 3 eV 2 [3l|. The value of the lightest mass, mi, was chosen to be 0.002 eV for the purpose 
of our numerical study, however there is no lower limit on the value of the mass of the 
lightest neutrino in either the normal or inverted hierarchies. If the mass of the lightest 
neutrino is lowered below the neutrino background temperature, Tc v b, then the resonance 
for this lightest state is governed by the corresponding thermal energy of the background 
neutrinos (see footnote 2 below). The particular choice for the neutrino masses in Eq. [21] 
results in the following features: 

1. When a scalar is produced, it decays predominantly to the m 3 mass eigenstate. Eq. 
[H gives P 3 « 0.967, P 2 « 0.031, and P y « 0.002. 

2. Because there is one scalar with mass Mq, Eq. [10] implies that the ratios of the 
neutrino masses govern the resonance energy positions, so that EJ^ S = 2/9 x E^ es and 
E* es = 1/25 x £f es . 3 

If we choose the lightest neutrino mass eigenstate to have the resonance at 15 MeV, to 
illustrate the effect, then the corresponding scalar mass is Mq 245 eV, and point 2 above 
implies that the other two resonances are E 2 es ~ 3.8 MeV and E^ es m 0.60 MeV. These 
two resonance energies are both well below experimental detection thresholds. Additionally, 
because the modified electron antineutrino flux is composed only of F Pl and Fp 2 the overall 
effect on F Pe is a depletion since decays of the scalar primarily contribute to F Pa . The 
resulting flux of electron antineutrinos can be seen in Fig. [1] Folding the flux with the 
cross-section for electron antineutrinos on protons, Eq. [20] gives the spectrum in Fig. [2] 

Relative to the electron antineutrino flux without interactions, the case with interactions 
has large depletion because of the dominant decay into the m% mass eigenstate. If we 
choose a heavier scalar so that the m2 mass eigenstate goes through the resonance at 15 
MeV instead of the mi mass eigenstate (i.e., Mq ~ 490 eV), then we get the same feature, 
but the depletion is smaller. This is because F Dl is multiplied by cos 2 9 12 ~ 0.70, while 
F 92 is multiplied by sin 2 12 ~ 0.30, so F 92 is the smaller component of the final electron 



3 The value of E^ es for the lightest neutrino could be in the range s/IEJTcvB % E^ es < \/2E v m\ 1 where 
the lower limit corresponds to the transition to the rclativistic case and Tc„b i* the background neutrinos 
temperature. Since TcvB ~ 2 x 10 -4 which is not far from mi (given that the effect goes like the square 
root of the mass in that range) our results will only be slightly modified when this is taken into account. 
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E v (MeV) 



5 10 15 20 25 30 35 40 

FIG. 1: The SRN electron antineutrino flux without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Major ana particles and for the normal mass 
hierarchy. 

a x F (events yr" 1 MeV" 1 ) 
17.51 




FIG. 2: The event rates for HyperKamiokande without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Majorana particles and for the normal mass 
hierarchy. 

antineutrino flux. If the scalar were even heavier (i.e., Mq ~ 1.2 keV) so that the m 3 mass 
eigenstate goes through the resonance at 15 MeV, then there is no depletion in the electron 
antineutrino flux since both the mi and rri2 eigenstates would have resonances at very high 
energies. As a result, neither of the corresponding fluxes would be visibly modified. 
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The mass of the lightest neutrino can also be lowered without bound (although as 
mentioned above we do not consider the cases where the mass of the lightest neutrino 
is below the mass of the cosmic background neutrino temperature). For a fixed value of Mq, 
as the mass of the lightest state is lowered, the position of the resonance moves to higher 
energies until it is in a region where the flux of the SRN is too small for the feature to be 
experimentally observable. 

To illustrate the effect on the electron neutrino spectrum, relevant to argon detectors such 
as Icarus and a future 100 kton argon detector lf|, we show the event rates in Fig. [3] for 
the normal hierarchy case. A next generation liquid argon detector with a size of 100 ktons 
could measure a significant number of electron neutrinos over just 5 years [if]]. At neutrino 
energies lower than about 19 MeV the solar neutrino flux dominates the SRN flux, and at 
energies greater than about 40 MeV the atmospheric neutrino flux begins to dominate. In 
Fig. [3j to obtain event rates, we have folded the SRN flux with the cross-section for electron 
neutrinos to interact with argon [if]]. We have used a resonance energy of 25 MeV for the 
lightest neutrino mass state (taken to be 0.001 eV so that Mq 225 eV). In this case we 
find significant reduction in the integrated event rate over the region in which the SRN flux 
is dominant. 




FIG. 3: The SRN electron neutrino flux for a normal hierarchy folded with the cross-section for 
electron neutrinos to interact with argon in a 100 kton detector running for 5 years. The red (solid) 
curve is for no interactions and the blue (dashed) curve is with interactions. 
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B. Inverted Neutrino Mass Hierarchy 



We now consider an example of the inverted mass hierarchy, characterized by mi 
m 2 >> 7713, with neutrino masses chosen to be 



m x = 0.05 eV, 
m 2 « 0.05 eV, 
m 3 = 0.008 eV. 



(22) 



311 ] for Am^i, the Vi-v 2 mass splitting. 



This reflects the best value of 7.92 x 10~ 5 eV 2 
Whenever a scalar is produced it dominantly decays into these two heavy eigenstates with 
equal probabilities. Because F Dl and Fp 2 are the contributing components to the electron 
antineutrino flux, in this scenario there are both regions of depletion but also regions of 
overall enhancement of the flux due to rescattering. We choose the case where the m\ and 
m 2 mass eigenstates go through resonance at 15 MeV, giving Mq « 1.2 keV, and show the 
results in Fig. H] (and Fig. [5] for weighting with cross-section). 



F (cm" 2 s" 1 MeV _1 ) 
3 




5 10 15 20 25 30 35 40 



E v (MeV) 



FIG. 4: The SRN electron antineutrino flux without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Major ana particles and for the inverted 
mass hierarchy. 

We find that the enhancement is large because all initial neutrino mass eigenstate fluxes 
produce scalars which add to the low energy m\ and m 2 eigenstate fluxes. The m\ and m 2 
eigenstate flux is depleted and is redistributed to lower energies. Once we fold the flux with 
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FIG. 5: The event rates for HyperKamiokande without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Majorana particles and for the inverted mass 
hierarchy. 

the cross-section, we see in Fig. [5] that in contrast to the case of the normal mass hierarchy, 
the peak at low energies is more pronounced. This is a result of the electron antineutrino 
flux being composed only of the m\ and m 2 neutrino mass eigenstates (since \U e ^\ 2 ~ 0), 
both of which get depleted by the resonance at the same energy since the two eigenstates 
are nearly mass degenerate. Therefore the electron antineutrino flux is almost completely 
depleted near the resonance cutoff. 

We apply similar analysis to the electron neutrino flux, of relevance to a liquid argon 
detector. In Fig. [6] we show event rates for lOOkton liquid argon detector. We use a 
resonance energy of 25 MeV for the two heavier mass eigenstates with mass 0.05 eV, which 
corresponds to Mq m 1580 eV. The integrated event rate for energies above the solar 
background are reduced relative to the no interaction case. 
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FIG. 6: The SRN electron neutrino flux for an inverted hierarchy folded with the cross-section 
for electron neutrinos to interact with argon in a 100 kton detector running for 5 years. The red 
(solid) curve is for no interactions and the blue (dashed) curve is with interactions. 

C. Quasi-Degenerate Neutrino Masses 

There still remains the possibility that the neutrino mass eigenstates are quasi-degenerate. 
For example, a mass hierarchy structure with 



(23) 



satisfies the requirements for the two independent mass splittings as well as cosmological 



mi r 


s 0.06 eV, 


m 2 f= 


0.06 eV, 


m 3 f 


a 0.08 eV 



constraints on the sum of the neutrino masses [32J. We consider the case where the two 
eigenstates with mass 0.06 eV have the same resonance energy as before, which would 
correspond to Mq = 1340 eV. The third mass eigenstate then has a resonance energy at 
approximately 11 MeV, however there is no corresponding dip since only the mi and ni2 
mass eigenstates contribute to the final v e flux. The result for the flux can be seen in Fig. 
[7] and for the event rate in Fig. [SJ 

It is clear by comparing Fig. H] and Fig. [7J that the case of the inverted mass hierarchy 
is nearly indistinguishable from the case of quasi-degenerate neutrino masses. However the 
quasi-degenerate case is distinguishable from the specific case of the normal mass hierarchy 
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E v (MeV) 



FIG. 7: The SRN electron antineutrino flux without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Major ana particles and for quasi-degenerate 
neutrino masses. 



a x F (events yr 1 MeV 1 ) 
17.51 




FIG. 8: The event rates for HyperKamiokande without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Major ana particles and for quasi-degenerate 
neutrino masses. 

when one neutrino mass is much lighter than the other two masses. 
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D. Multiple Depletion Dips 

There is the possibility within the normal hierarchy that the mi and m 2 mass eigenstates 
are nearly, but not exactly, degenerate, and also still much lighter than the m 3 mass 
eigenstate. One example of the possible values for the masses in such a scenario is 

mi pa 0.01 eV, 

m 2 ~ 0.013 eV, (24) 
m 3 ps 0.05 eV. 

If we choose the resonance of the m 2 eigenstate to be at 15 MeV, then the mi mass eigenstate 
goes through resonance at Ef es 20 MeV. This corresponds to M G ps 630 eV. 

As can be seen in Fig. [9] and Fig. [101 this leads to two depletion dips in the final 
electron antineutrino spectrum and three peaks. There is always a signal corresponding to 
each neutrino mass eigenstate interacting with the Goldstone. The presence of two distinct 
depletion dips in an experimentally interesting region, however, is sensitive to the ratio of 
the masses of two of the neutrino mass eigenstates, in this case mi and m 2 . For example, in 
Section IIII Al there is only one depletion dip because the ratio of the masses, and therefore 
the ratio of the resonance energies, for the mi and m 2 states is 4.5, so that with mi resonance 
at 15 MeV the m 2 resonance is at 3.8 MeV, outside of the observable region. 

Experimental observation of the energy position of these dips could determine the ratio 
of the mi and m 2 masses, which together with the measured value of Am^ and Am| 2 allows 
one to determine the neutrino masses. This is a remarkable possibility since it is extremely 
hard to experimentally determine the exact values of the neutrino masses, especially the 
mass of the lightest state. 

E. Dirac vs. Majorana Neutrinos 

If neutrinos are Majorana particles, then each scalar decay produces a or vrVr for 
each mass eigenstate. If the neutrinos are Dirac particles then the scalar can decay to vN 
or to Nv, where N and N are the extra neutrino fields added for the case of neutrinos being 
Dirac particles in the late-time neutrino mass generation models. Then only half of the 
decays of the scalar produce an antineutrino that will be seen in the detector. Therefore for 
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FIG. 9: The SRN electron antineutrino flux without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Majorana particles, for the normal mass 
hierarchy and where two neutrinos have distinct resonance features in the experimentally observable 
region. 

the case of neutrinos being Dirac particles there is an overall factor of 1/2 multiplying the 
last term of Eq. [T7] relative to the case of neutrinos being Majorana particles. 

If the neutrinos are arranged in a normal mass hierarchy as in Section UlI Al then the ability 
to distinguish between the neutrinos being Dirac or Majorana particles is confounded by the 
small amount of scalar decays into the m\ and m 2 eigenstates. However, in the case of 
the inverted mass hierarchy of Section IIII Bl there can be a visible difference in the electron 
antineutrino flux if the neutrinos are Dirac or Majorana particles as seen by comparing the 
Majorana particle case of Fig. H]with the Dirac particle case of Fig. [TTJ 

The resonance energy of the m\ and m 2 mass eigenstates have been set to 15 MeV in 
this case, exactly the same as for the Majorana particle case considered in Section IIII Bl 
Because of the extra factor of 1 /2 the overall scale of the enhancement is much smaller than 
in the case of neutrinos being Majorana particles. 



F. Non-adiabatic case (sin 



2 013 < 10- 5 ) 



If supernova neutrino flavor evolution is non-adiabatic, then the flux of electron 
antineutrinos that leaves a supernova is independent of the neutrino mass hierarchy (see 
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FIG. 10: The event rates for HyperKamiokande without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Majorana particles, for the normal mass 
hierarchy and where two neutrinos have distinct resonance features in the experimentally observable 
region. 
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FIG. 11: The SRN electron antineutrino flux without interactions (red (solid) curve) and with 
interactions (blue (dashed) curve) when neutrinos are Dirac particles and for the inverted mass 
hierarchy. 

Eq. [S] and Eq. ^ . If the supernova neutrinos interact with the cosmic background neutrinos 
via new light scalars then the flux observed at earth will be different for the normal mass 
hierarchy and the inverted mass hierarchy. However, this difference is not detectable at 
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the present or near-future neutrino experiments because it relies on the ability to observe 
neutrinos in the SRN flux at low energies where there is a large reactor background. 

The difference between the two neutrino mass hierarchies is present, even in the case of 
small sin 2 6 l i3, because for a normal mass hierarchy, the heavy neutrino mass eigenstate is 777.3 
which does not contribute to the v E flux. All of the scalars produced through the neutrino- 
neutrino interactions will dominantly decay into this heaviest neutrino mass eigenstate, and 
the final flux will have an overall depletion relative to the SRN flux without interactions. 
However, for the inverted hierarchy, the mi and m 2 mass eigenstates are the heavy states, 
while the mass eigenstate is the light state. The scalars produced through the neutrino- 
neutrino interactions will dominantly decay into the two heavy states, leading to a low 
energy enhancement as well as the higher energy dip. 

We show in Fig. [12] the SRN flux for the normal mass hierarchy (red/dotted curve) 
with m x = 0.001 eV, m 2 = 0.008 eV, m 3 = 0.05 eV, and M G = 173 eV, the 
inverted mass hierarchy (blue/dashed curve) with ni\ = 0.05 eV, rri2 = 0.05 eV, rris = 
0.008 eV, and Mq = 1225 eV, and the SRN flux without interactions (black/solid curve). 
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FIG. 12: The SRN electron antineutrino flux without interactions (black solid curve), with 
interactions and normal mass hierarchy (red dotted curve) , and with interactions and inverted mass 
hierarchy (blue dashed curve), when neutrinos are Majorana particles and for sin 2 6\3 < 10 -5 . 
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G. Signal Detection 

Here we show an example of an inverted mass hierarchy. The resonant energy for the 
mi and neutrino mass eigenstates is taken to be 16 MeV, which for m\ ~ m>2 ~ 0.05 eV 
gives Mg ~ 1265 eV. We present both the expected SRN flux as well as the SRN flux folded 
with the cross-section given in Eq. [20] in Fig. [13] and Fig. [T4] 




FIG. 13: The SRN electron antineutrino flux without interactions (red solid curve), and with 
interactions and inverted mass hierarchy (blue dashed curve) when neutrinos are Majorana 
particles. 



Comparing the blue (dashed) curve to the red (solid) curve in Fig. [THwe see that there is 
a depletion of the SRN event rate above approximately 8 MeV if the SRNs interact via the 
light scalar at resonance with the cosmic background neutrinos (blue dashed curve). This 
depletion is present up to the location of the dip, which in this case is 16 MeV. The main 
source of background, assuming the addition of Gd to the water Cerenkov detector, are 
nuclear reactor electron antineutrinos, but this background is small above neutrino energies 
of about 12 MeV Q, Q]. This background is dependent on the location of the experiment 
and could be nearly absent. Clearly the detector which is not near nuclear reactors would 
have a better chance of seeing the signal for lower values of the resonance energy HQ. 

To demonstrate the significance of our signal we look at an energy of 15 MeV (where the 
effect is most significant and the reactor background is negligible). At this energy we expect 
approximately 11 events per year at HyperK from the SRN flux without new interactions 
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FIG. 14: The event rates for HyperKamiokande without interactions (red solid curve) and with 
interactions (blue dashed curve) when neutrinos are Majorana particles, for the inverted mass 
hierarchy. 

(red solid curve). However, if the neutrinos interact via the light scalar at resonance, then 
we instead predict approximately 2 events per year (blue dashed curve) at this energy bin 
(assuming a 2 MeV bin). The average fluctuation in the number of events expected in 5 
years with no interactions can be estimated as a ~ v55 ~ 7.5 events. If we take as our signal 
the number of events expected without interactions minus the number of events expected 
with interactions (i.e., the event deficit) over the 5 year period (this is 55 — 10 = 45 events), 
then in 5 years one expects approximately a 6cr (45/7.5) effect. A similar analysis can be 
performed for the previous cases discussed in this paper, but since the depletion in these 
cases is at lower energies one must pay careful attention to the reactor background. 

This analysis also requires a side band study, in order to determine the SRN flux in 
a region where interactions are ineffective (in our present example, this would be above 
16 MeV). This would provide the overall normalization necessary for establishing the 
background. Clearly, once the shape of the signal without resonance and the non-SRN 
background are known, a more sophisticated analysis (including a bin-by-bin fit to the 
shape of the curves) can be achieved, and may even provide enhanced significance. This, 
however, is beyond the scope of this work, whose aim we regard to be a discussion of the 
qualitative aspect of our new physics signal. 

For a liquid argon detector, if the resonance energy is above the cutoff for the solar 
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neutrino flux at about 19 MeV, the integrated number of SRN electron neutrino events 
would be visibly reduced in the presence of new interactions. This is shown for a normal 
mass hierarchy and an inverted mass hierarchy in Fig. [3] and Fig. [6] respectively. While our 
total number of events is a conservative estimate (we use a z evolution that is flat above 
z — 1, instead of stronger dependence 16] ), the depletion of events up to the resonance 
cutoff is a robust feature. We find that there is approximately a 25% reduction in the 
number of electron neutrino absorption events with new interactions compared to without 
new interactions for both the normal and inverted mass hierarchies for the energy window 
from 19 MeV to 40 MeV. 



IV. CONCLUSION 

The late-time neutrino mass generation models could be tested by detecting unique 
features of the SRN flux in both its electron antineutrino and neutrino components (for other 



tests of new physics that can be done with the SRN flux see 33J). To illustrate this new 
effect we have considered an Abelian U(l) model that generates neutrino masses at low scales. 
However, it is clear that the main features would still hold for a more complicated, non- 
Abelian model, although these models are already more constrained by BBN considerations. 
For example, one could still have observable dips in the SRN spectrum if the resonances are 
in a desirable energy window, but the couplings are no longer proportional to the neutrino 
masses, and so some predictive power is lost. However, one could hope to correlate the 
observations of the dip locations in the SRN spectrum with signals proposed to be present 
in the cosmic microwave background [?| for this case. We expect that the future generation 
water Cerenkov detectors enriched with gadolinium such as UNO, HyperKamiokande, or 
MEMPHYS would be able to detect a substantial number of SRN antineutrino events in 



a year 15j. Note that the threshold for this is on the order of 10 MeV and depends on 

nn 

the location of the detector, especially due to reactor backgrounds [13J, |14J. The effects of 
smearing due to the energy resolution of the water cerenkov-type detector needs to also be 
taken into account in a detailed analysis. For a Gaussian energy resolution function with 



width A (A/MeV ~ 0.6 y/E/MeV), the smearing is only at most a few MeV [13] in the 
energy domain we considered. This smearing is then always smaller than the width of the 
depletion features considered. The neutrino component of the SRN flux could be detected 
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by a large 100 kton liquid argon detector 16J. If there are neutrino- neutrino interactions 
through the light scalars present in these models, there is a possibility to distinguish between 
normal and inverted mass hierarchies and Dirac versus Majorana neutrinos, as well as to 
determine the absolute values of the neutrino masses. The ability to distinguish the neutrino 
mass hierarchy is independent of whether supernova neutrino flavor evolution is adiabatic 
or non-adiabatic. 

The qualitative features of the signal of new interations via light scalar, such as the 
depletion, enhancement at lower energies, and the possibility to distinguish between neutrino 
mass hierarchy, as well as the nature of the neutrino are independent of the theoretical 
model for the supernova neutrino energy spectrum, which predict slightly different shape 
and wider range of average energies for different neutrino (antineutrino) flavors than the 



KRJ model 
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21] . This is not surprising because the produced neutrino spectrum 



does not depend on the detailed shape and normalization of the initial supernova neutrino 



fluxes but rather on the coupling of the sea 



ar to the final neutrino mass eigenstates. As 



shown in section A (3) of the first paper ll], the presence of a deep dip is universal, and 



its position depends only on the masses of the scalar and the target neutrino, rather than 
any feature of the neutrino spectrum. For a normal neutrino mass hierarchy there is an 
overall depletion of the SRN flux, while for an inverted neutrino mass hierarchy there is an 
enhancement of the SRN flux at low energies and a region of depletion at higher energies. 
If a sterile neutrino with a much larger mass than the active neutrinos were also to couple 
to the new scalar, then independent of the details of the masses of the active neutrinos the 
effect would be almost complete depletion of the spectrum in some energy window since the 
scalar would decay predominantly to the massive sterile neutrino. 

All of these signals, and especially their observation, depend on the parameters of the 
model. In Fig. [15] we show constraints on the parameter space for which the SRN effects 
can be obtained in the y v — Mq plane. The signals proposed here are present in the SRN 
flux only if the couplings of the neutrino mass eigenstates to the scalar are larger than the 
condition given in Eq. [12] for a given value of Mq. This condition comes from requiring 
that the mean free path for absorption of a SRN neutrino on a cosmic background neutrino 
is much smaller than the Hubble scale 11|. It is a sufficient condition to guarantee the 



absorption of all three neutrino flavors. This lower bound on the coupling is represented 
by the diagonal blue (solid) line. If the resonance energy is below 12 MeV then there is a 
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large background from nuclear reactor antineutrinos [13I. \h\. To have a significant signal we 



take M G to be above J 2m u E^ in , where E^ in is approximately 15 MeV. These threshold 
values are represented by the three vertical red dashed lines which are calculated for values of 
m v = 0.001 eV, 0.008 eV, and 0.05 eV. If the mass of the scalar is larger than these values, 
then the signal would be above the reactor background. Similarly, the signal would not be 
observable if the mass of the scalar is large so that the heaviest neutrino mass eigenstate has 
a resonance energy in the region where the SRN flux is small. We also show the constraint 
imposed by BBN considerations, which is similar to the bound obtained from SN cooling 
and to the bound from the observation of undegraded SN1987A neutrino flux The 
SRN flux is also sensitive to the non-resonant process, for example 2v — > <p — > 2G — > 4u, but 
only in a very small region of the parameter space, above the horizonal black dashed line 



and below the horizontal red solid line 



111 ]. The area above the diagonal green dashed line 



corresponds to the BBN constraint for a non-abelian Majorana case. We note that there is 
still a large range of parameter space where the couplings are large enough to give SRN flux 
modification in an energy window that large neutrino detectors could directly probe. 

We have shown that the cosmic background neutrinos interacting with supernova relic 
neutrinos through exchange of the light scalar lead to significant modification of the SRN 
flux observed at earth. These signals would be detectable for a large region of parameter 
space, in some cases at a significance of more than 5a, and measurements of the presence 
of these effects are well within the reach of the next-generation water Cerenkov detectors 
enriched with gadolinium, or a large 100 kton liquid argon detector. Specifically we have 
shown that the changes induced in the flux by the exchange of the light scalars might 
allow one to distinguish between neutrinos being Majorana or Dirac particles, the type of 
neutrino mass hierarchy (normal or inverted or quasi-degenerate), and could also possibly 
determine the absolute values of the neutrino masses. An interesting feature is that the 
ability to distinguish neutrino mass hierarchy does not depend on the dynamics of the flavor 
evolution of neutrinos leaving the supernova (whether it is adiabatic or non-adiabatic), or 
on the specific shape and normalization of the initial supernova neutrino flux. Note that the 
hierarchy determination can be made by solely looking at the spectrum of supernova relic 
electron antineutrinos, without need to do the measurement of the flux of supernova relic 
electron neutrinos. In addition, the modification of the SRN flux in any of the proposed 
scenarios is a clear indication of the presence of the cosmic background neutrinos left over 
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from the era of Big Bang Nucleosynthesis. 
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